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Abstract 

The differential transformation method (DTM) enables the easy construction of a power-series 
solution to a nonlinear differential equation. The exponentiation operation has not been specif- 
ically addressed in the DTM literature, and constructing it iteratively is suboptimal. The recur- 
rence for exponentiating a power series by J.C.P. Miller provides a concise implementation of 
exponentiation by a positive integer for DTM. An equally-concise implementation of the expo- 
nential function is also provided. 
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1. Introduction 

Constructing power-series solutions to differential equations, especially those which do not 
^ . admit a closed-form solution, has long been an important, and widely-used, solution technique. 

00 ' Traditionally, computing power-series solutions required a fair amount of "boiler-plate" symbolic 

C^^ i manipulation, especially in the setup of the power-matching phase. The differential transforma- 

tion method (DTM) enables the easy construction of a power-series solution by specifying a 
conversion between the differential equation and a recurrence relation for the power-series coef- 
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The table in the current literature which specifies the translation between the terms of the 
differential equation and the recurrence relation has a striking omission: it contains no expo- 
nentiation operation. Exponentiation by a positive integer can be constructed iteratively using 
the table entry for multiplication (i.e. multiplying the function with itself n times), but such a 
\^ . construction is suboptimal because it leads to n - 1 nested sums. Using a recurrence for expo- 

V^ ' nentiating a power series by J.C.P. Miller, a table entry for positive-integer exponentiation can 

be provided which introduces only a single sum. A single-sum recurrence for the exponential 
function of a power series can be similarly constructed. 



Email address: hal . f inkelSyale . edu (Hal Finkel) 



Preprint submitted to Elsevier July 14, 2010 



2. Preliminaries 



2.1. Differential transformation 

The differential transform of the function w(x), called W(k), is defined as Uj] 
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where jj is the k derivative with respect to x. The inverse transformation is 

w(x) = ^ W(k)x^ . 

k=0 



(2) 



Combining Equations[T]and|2]it is clear that the differential transform is derived from the Taylor- 
series representation of w{x) about x = 
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dx^ 



w{x) 
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Table [T] shows the traditional list of operations and their differential transform^ 



Original Function 


Transformed Function 


w{x) - y{x) + z(x) 

w{x) - Ay{x), A a constant 

w{x) = §^y{x) 

w{x) - y{x)z{x) 

wix) = JC™ 


W(k) = Y(k) ± Z{k) 
W{k) = AY{k) 

W{k)^Y.UY(l)Z{k-T) 

W(k) = 6k,m, 6 is the Kronecker delta 



(3) 



Table 1 : Traditional operations under DTM 

2.2. J.C.P. Miller's recurrence 

In Knuth's classic book on seminumerical algorithms \^, he provides a recurrence relation 
attributed to J.C.P. Miller for the coefficients of a power series raised to some integer power 
Unfortunately, the recurrence is not well known \^. Given some formal power series 
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k=a 



where A^ could be oo and aq ^ 0, then 
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' Recent literature has also included tables with entries for additional elementary functions (e.g. sin(A;), cos(a:), e') and 
integral relations \^ 



where the {q} are given by the recurrence relation j 



C() = a; 







Q = I^ llj=i[(m + 1); - kjOjCk^j 
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A simple proof by Zeilberger is as follows Ja]: Ck is the coefficient of x in the Laurent series 
expansion of w(x)"^/x''. For any Laurent series, fix), the coefficient of x" in x-^f{x) is zero. So: 
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[x"] I — A:(flo + aix + . . . + a^x^) \ 

+(m + l)(ai + 2a2X + ... + Na^x^-^ ^) 
= [x'] {-k[ao^+a,^ + ... + a^^] 

+(m+l)[a,^+2a2^ + ...+Na!,^]) 
= -k[aoCk + aiQ-i + . . . + flA?Q-A?] 

+(m + l)[aiQ_i + 2fl2Q-2 + . • . + NuNCk-N] 

From which the recurrence follows. The standard notation ISfl that [x"]/(x) is the coefficient of 
x" in fix) has been used. 

3. Exponentiation in tlie DTM 

Combining the recurrence relation with the DTM formaUsm is straightforward. Note that 
W{k) is the the k Taylor-series coefficient of w(x), and this yields the table entry in Table ID 
Care must be taken, however, to insure that if F(0) is zero then the replacement ^(x) -^ xy(x) 
where y{x) - — is made. 



Original Function 


Transformed Function 


w(x) = yixf, F(0) ^ 
w(x) = e><^> 


"^ik) = j^ i:%,{{m+l)j-k]Y{j)W{k-j) 

W(0) = 7(0)"' 

W{k) = \ 2 ■-! jY{i)W{k - j) 

W(0) = e^C' 



Table 2: Exponentiation operations under DTM. n is a positive integer 

Using this formulation greatly simplifies the expressions resulting from exponentiating a 
function compared to the previously-available method yl]. Specifically, for w(x) = ^(x)'" the 
representation constructed iteratively from the multiplication rule was 



oo k /;,„ 
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^(^^ = Z Z Z ■ ■ ■ Z Z ^^^i)^(^2 - ^i) ■ ■ ■ F(^m-i - k„^2)Yik - fe„,_i) . (8) 
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Numerically, there is a significant advantage to using the recurrence in Table |2] which requires 
O(A^) computations, compared to the iterative method (Equation[8]l, which requires 0(mNlogN) 
computations. 



4. The exponential function in DTM 

Using Zeilberger's method, it is possible to derive a recurrence for e"'*^^: 

= [A {-k^ + («i + 2«2X + . . . + Na^x^-' 5^) ^^^ 

= U°] {-k^ + [fli §^ + 2fli g^ + . . . + A^flAT^]) 

= -^Q + [fliQ_i + 2fl2Q-2 + . . . + Na^Ck-N] 

The resuhing recurrence is also listed in Table |2] Again, it is much simpler than the iterative 
exponentiation construction in combination with the Taylor-series expansion of the exponential 
function ^ 

oo ^ iM k k,„-i k^ ^2 

m=o'^' Jl=0yfc„,_i=0/t„,_2=0 /(:2=0y(:,=() 

5. An example: The One-Dimensional Planar Bratu Problem 

To exemplify the formulation presented here, the simplified recurrence will be applied to the 
one-dimensional planar Bratu problem which had previously been studied using DTM by Hassan 
and Ertiirk [2]. The differential equation is 



u" + Ae" =0 < jc < 1 

m(0) = 0,m(1) = 



(11) 



and the previously-derived DTM solution is 



(k+l)(k+2) 
U(ki)Uik2 -h)--- f/(fc„,_i - k,„^2)U(k - k„,^i)) 



for which the coefficients {U{k)} are computed starting from 

t/(0) = 

c/(i) = r 



(13) 



where y is a constant fixed by the x = 1 boundary condition 

oo 

J]u(k)^0. (14) 

This can be compared to the analytic solution yj] 

cosh((x-i)f) 



u(x) = -2 In 
where solves 



(15) 



cosh(f) 
^/2Acosh(-\ . (16) 



Depending on the value of A, Equation [16] has zero, one or two solution(s). See the above- 
referenced paper by Hassan and Ertiirk, and references cited therein, for further details. 
Applying the recurrence for e"^^' provided here yields the following DTM solution 

W(k)^lY.%iJU(j)W(k-j) (17) 

W(Q) = e^C' 

which can be simphfied by writing W(k) = - »+')(^+2) Uik + 2) to 
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U(k + 2) = «iTiW2) 27=1 M - J + 1)(^ - ; + 2)Uij)U(k -j + 2) k>l 



k(k+\)(k+T, ^MJ"-"^ 7 -^ ^yw j-r^,^yj,^y^ J -r ^j ^ i. i 

U{2) = -^^V{Q) ^^»^ 



The first few values of the recurrence are 

U{Q) = 

t/(2) = -f ^^^^ 

f/(3) = -rs 

and these are the same regardless of whether Equation [12] or Equation [TSl is used. 

6. Conclusion 

The rule derived from J.C.P. Miller's recurrence, and a similar recurrence derived for the 
exponential function, given in Table|2] are far simpler than the representations previously used in 
the literature (Equations [8] and [TOll. It should prove easier to apply DTM, and any power-series- 
solution technique, to nonlinear differential equations using the recurrence relations presented 
here. 
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